Exponential Growth of Particle Number far from the Parametric Resonance Regime 
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Parametric resonance has received a considerable amount of interest as a good mathematical model 
to describe the initial stages of the reheating phase (matter creation) in inflationary cosmology. It 
is also known that exponential particle creation can occur in situations which do not fall in the 
parametric resonance regime characterized by oscillations of the inflaton field about its minimum. 
Here we present a new analytical approach to exponential particle production which can occur 
when the inflaton is far from the minimum of its potential. Crucial for this effect is a term in the 
equation of motion which acts like a negative mass square term, as occurs for tachyonic preheating 
and negative coupling particle production. Our techniques apply in models with a strong coupling 
between matter fields \ an d the inflaton <f>, or in some models in which the inflaton has a large 
amplitude of oscillation. Note that our analysis yields results which are quite model dependent. 
Exponential growth occurs in a model with interaction Lagrangian — gA4 p i(f>x 2 - However, for the 
interaction Lagrangian —g 2 4> 2 X 2 i our formalism shows that in the large coupling limit there can only 
be exponential particle production when <f> crosses 0. 



I. INTRODUCTION 

Inflationary Cosmology is currently the most widely 
accepted model of the very early Universe. Inflation is 
expected to end in a period in which the inflaton field <p is 
evolving rapidly and, finally, oscillating about the ground 
state of its potential. During this phase - commonly 
called the reheating phase -, energy is transferred from 
the inflaton field to ordinary matter, here modelled by 
another scalar field \, via processes which depend on the 
coupling between <fi and \ m t ne interaction Lagrangian. 
The details of this energy transfer are important for sev- 
eral key issues in cosmology, among them baryogenesis, 
the relic gravitino abundance, and the possible produc- 
tion of topological defects. For example, in order for the 
GUT-scale baryogenesis mechanism to be effective, the 
temperature of matter after reheating has to be compa- 
rable or higher than the GUT scale T = 10 16 GeU (for 
recent reviews see e.g. [Q ) . If the reheating process is slow 
(i.e. the decay rate of the inflaton is small compared to 
the Hubble expansion rate at the end of inflation), as is 
predicted using a perturbative analysis of inflaton decay 
P,|l, then the temperature after reheating is predicted 
to be much lower than the GUT scale. 

Recently, however, it was realized that in many 
models of inflation the energy transfer takes place much 
more rapidly, based on the phenomenon of parametric 



resonance. The basic idea is the following: while the in- 
flaton field is oscillating about its minimum, it induces a 
periodically varying time dependence in the mass of the 
matter field, which in turn leads to exponential growth of 
the matter fields which corresponds to an explosive cre- 
ation of particles. The theory was put on a firm mathe- 
matical basis in (these analyses included the expan- 
sion of the Universe), and a comprehensive and detailed 
study was given in |j. According to the new picture, the 
whole reheating process is divided into two phases. The 
first one is the phase of parametric resonance, called pre- 
heating Q. This is an out-of-cquilibrium process during 
which the growth of the number of particles is exponen- 
tial but only occurs for momenta in certain resonance 
bands. In the second phase, interactions between the 
matter particles and back reaction effects become crucial 
and lead to thermalization of matter Q . Although para- 
metric resonance is an attractive mechanism for particle 
production, there are other possibilities to obtain expo- 
nential particle production, for example instant preheat- 
ing 0, chaotic preheating |To| and tachyonic preheating 
|1 l| , |L2| . In instant preheating, particles are produced ex- 
ponentially when the inflaton field (f> crosses the minimum 
value and the interaction potential is negative (see also 
p3[ for a discussion of negative coupling instability) , in 
chaotic preheating the exponential increase in the per- 
turbations is driven by the chaotic dynamics of the back- 
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ground, and tachyonic preheating is particle production 
driven by a negative effective mass square term which 
occurs during spinodal decomposition in models of new 
and hybrid inflation |3 . 

In this paper we present a new mathematical approach 
to exponential particle production at the end of inflation 
but far from the parametric resonance regime. Our ef- 
fect is generated by a mass square term for the matter 
fields which is negative for certain periods in the evolu- 
tion of the inflaton field (in this respect the mechanism 
is exactly what happens in tachyonic jll],[[2| and neg- 
ative coupling 13 preheating). However, our analytical 
method applies also for quite general nonlinear correction 
terms in the equations of motion. Note that we are not 
assuming periodic variation of the external field. Typi- 
cally, our mechanism operates for large values of \(f>\ but 
only for particular forms of the interaction Lagrangian (in 
particular we need a large effective coupling constant). 
We derive our effect making use of some results on the 
asymptotics of solutions of second order ordinary differ- 
ential equations Jl5| . 

In next section we summarize the relevant mathemat- 
ical results quoted from J15|. In Section III we apply 



them to differential equations which describe the evolu- 
tion of the matter fields during reheating in the presence 
of a dynamical inflaton field. Finally, we briefly discuss 
the results and point out analogies and differences with 
previous work. 



II. MATHEMATICAL APPROACH 

The results below are based on the asymptotic analysis 
of second order differential equations as described e.g. in 
reference [jL5[ . 

We consider the ordinary differential equation 



V - p(t, o)y = , 



(1) 



where p(t, a) is a real function for values of t in some 
interval /, and for values of a in an interval E. The 
derivatives are with respect to t. Both t and a are real 
variables, and in our applications to cosmology t will be 
time, whereas a is some external parameter. In this sec- 
tion we recall some results concerning the behavior of the 
solutions yi(t) and 2/2 (0 of (|l|) for values of t in I when 
a — > +00. 

We assume that in the limit a — > +00, the function 
p(t, a) can be expanded in the following way: 



p{t, a) 



n=0 



(2) 



where r is a positive rational number. We work in subin- 
tervals /' C I where po(t) does not vanish, but could have 
at most double poles in I', 

With these hypotheses, by a theorem in reference Jl5| , 
it can be shown that: 



a) In subintervals I' a I where po(t) > 0, and when 
— ► 00 : 



yi(t) = {powr 1/4 ex P |y K /2 ( Po (i')) 1/2 + 

l pi (t')(po(0)- 1/2 ]^'}{l + o(^} (3) 



and 



= {PoWr 1/4 ex P |- J [a r / 2 (p (t'))V 2 + 
i pi ^)bo(i'))- 1/2 ]rfi'}{l + 0^)| (4) 

b) In subintervals I' C I where po(t) < 0, and when 
(7 — ► 00 : 



y 1 (t) = c{-p (t)}- 1 / i S m\ I [W 2 (-p (t')) 1/2 + 
1 



Pi(t')H>o(Or 1/2 + a. 



]^}{i + o(i)} 



(5) 



and; 



y 2 (t) = c{-p (t)}- 1/4 cos| J [W 2 (-p (O) 1/2 + 

i pi (0(-rf)r 1/2 + a]dt' }{l + ( 6 ) 

Here, the phase angle a is a real constant, and c is a 
possibly complex constant. 

Thus, if po(t) > on J', in the asymptotic approxima- 
tion y(t) increases (or decreases) exponentially. However, 
iipo(t) < 0, we obtain oscillating solutions. This theorem 
generalizes the trivial dynamics of a harmonic oscillator 
with positive or negative square mass to the case where 
the coefficient representing the square mass is not con- 
stant in time (but does not cross zero, either). Thus, in 
the adiabatic limit in which the mass varies only slowly, 
the results are easy to understand from the point of view 
of physics. These results will be used in the next section 
in order to study exponential growth of particle number 
at the end of inflationary phase of the evolution of the 
Universe. 

III. APPLICATIONS 

In this section we use the above theorem to get some 
information on particle creation induced by the inflaton 
field. In order to be as general as possible, we do not 
prescribe any specific equation of motion for the inflaton 
field, but only consider some general properties needed, 
which are shared by several models of chaotic inflation. 
As indicated, in the following 4> denotes the inflaton, and 
X a scalar matter field. As in most studies of preheat- 
ing, we will neglect the bare mass and self-interactions 
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of the x field. We will also neglect the expansion of the 
Universe. As long as the time scale for exponential in- 
crease obtained under this approximation is smaller than 
the Hubble expansion time, we expect the approxima- 
tion to be a good one. As interaction Lagrangian we 
take Li n t — ~5/(0)x 2 j where the function / contains a 
large parameter which plays the role of a in the previous 
section. This parameter can be a coupling constant - it 
can also be the amplitude of the inflaton during a period 
of oscillation. 

The equation of motion of the k-th Fourier mode of 
the scalar matter field x is given by : 

Xk - (-fc 2 - f(4>)) X h = o . (7) 

Comparing this with equation (|l|) we see that: 

p{4>) = -fc 2 - f(<p) . (8) 

As mentioned above, we can write f(<fi) — ah(<fi), 
where a is a large parameter. Then we can write the 
above equation as 

p(0) = -fc 2 - ah{4>) . (9) 

Comparing with the expansion (||) we get 

po(t) = -ah{<f>{t)) ; 

p 2 (t) = -fc 2 ; (10) 

p n (t) = 0; n^0,2 

Now, from the key theorem in section [n] it follows 
that if the parameter a is large, then in time intervals 
where h(<f>(t)) is positive the solutions exhibit exponential 
growth, whereas in those during which h(cj)(t)) is negative 
we get bounded oscillating solutions. 

In order to connect the above results with the physics 
of known models we now specialize to two cases: 

A. Case /(</>) = <? V 

In this case we can write p(4>) as 

p{4>) = -fc 2 -.g 2 2 <I> 2 (i) (11) 

where $ = and 0o is the amplitude of the infla- 
ton field at the beginning of the oscillatory phase during 
which parametric resonance occurs. Taking a — g 2 <^\ 
and r = 1 in the expansion (^|) we have: 

po = -$ 2 (t) ; 

P2 = -fc 2 ; (12) 
p n = ; n ^ 0, 2 

Since in this case Po(t) < for all t, the solutions Xk 
are always asymptotic oscillating when a is big enough. 
Thus, in this model no particle production occurs dur- 
ing time intervals when <f> ^ 0. This result sheds new 
light on the observations QJ^] that particle production is 
concentrated at times when <j) = 0. 



B. Case f(<j>) = gM p i<}> 

In this case: 

p{<t>) = -fc 2 - gMjaftt) . (13) 

We assume that the potential can be well approximated 
as quadratic near the origin, and as linear for large val- 
ues of \4>\ (|0| > Mp{). The second assumption is made 
in order to obtain slow- rolling of <j> for large values of \<j>\. 
In an improved analysis the slow-rolling phase would be 
generated by the expansion of the Universe and the re- 
sulting Hubble damping term in the scalar field equation 
of motion. Since 4>o is of the order M p i, the condition to 
apply the theorem of Section (||) (to obtain large a) is 
fc 2 ,m 2 <C gMpifio, where m is the inflaton mass. Thus, 
in particular, 

' » te) 2 

For a potential of the form V(</)) = rn 2 <j) 2 + A</> 4 the 
COBE constraints on the fluctuations generated during 
inflation give m/M pl < 1CT 6 and A < 1CT 12 . Thus, @ 
becomes g ^> 10~ 12 . On the other hand, by naturalness, 
g < vA. Therefore, the large coupling asymptotic analy- 
sis is valid in a large interesting region of values of g from 
10" 12 to 1CT 6 . 
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FIG. 1. Sketch of a chaotic inflation potential V(4>). 
[$| > 1 is the slow rolling phase in which a linear approx- 
imation to V{4>) is used. The region |$| < 1 is the region of 
oscillation with a quadratic approximation. 

Taking a = gM v i<j>o and r = 1 in the expansion (j^), 
and using the same variable $ introduced in the previous 
subsection, we have from (|l3|): 

Po (t) = -*(*) 

P2 = -fc 2 (15) 
p n = ; n ^ 0, 2 . 

It follows immediately that for $ < we have po(t) > 0, 
and by the theorem presented in Section || we expect 
exponential particle production. 
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In the following we assume that <fi starts out at large 
negative values in the slow-rolling phase. Once it reaches 
|$| ~ 1, 4> will proceed to oscillate about the minimum 
of the potential. Note that exponential increase in x is 
possible in both regions, as long as (j> < 0. In the following 
we analyse the two phases separately. In both phases, the 
equation of motion is 



d 2 ^ dV{(f) 



. 



(16) 



1. Slow rolling region 

In the slow rolling phase we take a linear approxima- 
tion to the potential, i.e., V{4>) — A<fi + B, where A and 
B are constants. In the region <p < ~4>q 1 the constant 
A is negative, and B is chosen such that the potential is 
continuous at the transition points |$| = 1. 

Solving Equation (Jlq) in this approximation gives: 



Ct + D 



(17) 



where C and D < are constants. Thus, 



p{(f) ) = -k 2 -g^4t 2 - g 2 Ct-g 2 D 



(18) 



Taking a = g 2 \A\/2 and r = 1 in the expansion 
we obtain from (Uq): 



Po(t) = 
P2 = -k 2 ; 
Pn = 0;n^0,2 



{t 2 + Ct + D); 



(19) 



where C = 2C/\A\ and D = 2D/\A\. 

There are two roots t- and t+ of the equation po = 0: 



t± 



-C± 



C 2 -AD 



(20) 



Since D < 0, the roots are real. Thus, there is a time in- 
terval t + ] during which the solution for Xk is asymp- 
totically exponential. In fact it is easy to see that C < 
and At = t+ — t_ > \C\ = —C. From this inequality we 
conclude that the smaller A or the bigger C, the bigger is 
the time interval during which the solution has exponen- 
tial growth. The case of A small suggests that chaotic 
inflation, with its long slow roll-over section of the poten- 
tial, is the most effective for particle production, when 
the system is far from the parametric resonance regime.. 

Note that in the slow rolling region with (j> > <j>o , po (t) 
is negative and, as expected, we obtain asymptotically 
oscillating solutions. 



2. Oscillatory region 

For |$| < 1 we take a quadratic approximation, that 
is, V(4>) = A'(f) 2 , where A' > 0. Solving the differential 
equation (|l^ ) for 4> we get 4>{t) = <j> cos(wt), where 4>q 
is the amplitude of oscillation and w = \[2A. Thus we 
have: 



p = —k 2 — gM p i4>o cos(wt) 



(21) 



Taking again a = gM p i(f>o and r = 1 in the expansion 
j|) we obtain from ( pl| ) 



Po = — cos(wt) ; 

P2 = -k 2 ; 

p n = ; n £ 0, 2 



(22) 
(23) 



For (j — iW < wt < (j + iW with j even, we have 
Po < 0. So if a is large, which can occur either if the 
coupling constant or the field amplitude are large, the 
theorem of Section || give us an oscillating asymptotic 
solution: 



Xk(t) = 



^/cos(utf) 



: COS 



[(7 1/2 y/cOs(wt') + a]dt' 



(24) 



where a and c are constants. 

However, for time intervals with (j — ij7r < wt < 
(j + |)7T with j odd, we get po > 0. Then, for large a, 
the asymptotic solution reads 



1 



Xk 



yf— COs(wt) 



t 1/2 J* y/- cos(wt')dt' 



(25) 



and in this case we get exponential production of matter. 

Combining both cases, we get the time evolution of the 
particle number density rife drawn in Figure (||). Here, nk 
is the number density of particles produced in the k-th 
mode obtained using the standard formula (see e.g. 0) 



n fc = w k 



\Xk\ 



\Xk\ 



1 

2 ' 



(26) 



where Wk is the effective frequency of Xk- Substituting 
( p4| ) and ( p5[ ) in ( p6| ) within their respective intervals of 
validity, we obtain (as a function of time) a kind of pulsed 
particle creation (see Fig. 2). Note that our results in 
this region are similar to what was previously obtained 
in the context of tachyonic preheating |0| . 



IV. DISCUSSION AND CONCLUSION 

We have presented a new analytical approach to expo- 
nential particle production far from the parametric res- 
onance regime. Our analysis applies to situations when 
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In n k 




FIG. 2. Pulsed Particle Creation 

the either the absolute value of the inflaton field is large, 
or the coupling constants are very large. The results, 
however, are is quite model-specific. One obtains expo- 
nential particle production for the interaction Lagrangian 
L int = -gMpi^x 2 , but not for L int = -g 2 4> 2 x 2 - Crucial 
to have exponential particle production is the presence of 
a negative mass square term in the effective potential for 
the inflaton. Note that our method is robust in the sense 
that the results do not change qualitatively when adding 
extra terms in the equation of motion such as explicit 
mass terms and self interaction terms for the x field, as 
long as the additional terms are small in amplitude in 
the appropriate units. 

In the rolling phase we can have exponential particle 
production when the coupling constant or the slope of 
the potential are large. However, the length of the time 
interval where this exponential growth happens decreases 
when the slope A of the linear potential V(4>) = A<f> + 
B increases. Then, when A is large enough to use the 
asymptotic approximation, we do not have exponential 
matter creation for a long time. On the other hand, if A 
is small (yielding a large interval of exponential growth) , 
we only can use the asymptotic approximation provided 
the coupling constant g is big enough. In this case we 
can have particle creation for a long time interval even 
during the rolling period. Applied to the regime of fast 
rolling, we are in the situation of tachyonic preheating, 
and our results agree with those of py,F2[ ■ 
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